On thermal stability of topological qubit in Kitaev's 4D model 
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We analyse stability of the four-dimensional Kitaev model - a candidate for scalable quantum 
memory - in finite temperature within the weak coupling Markovian limit. It is shown that, below 
a critical temperature, certain topological qubit observables X and Z possess relaxation times 
exponentially long in the size of the system. Their construction involves polynomial in system's 
size algorithm which uses as an input the results of measurements performed on all individual spins. 
We also discuss the drawbacks of such candidate for quantum memory and mention the implications 
of the stability of qubit for statistical mechanics. 



INTRODUCTION 



While quantum computation oflFers algorithms which 
can outperform the classical ones, they are very fragile 
with respect to external disturbance. Therefore, along 
with the discoveries of fast algorithms, the question of 
how to protect quantum computation against decoher- 
ence was the subject of extensive studies. As a result the 
whole domain was created known as fault tolerant quan- 
tum computation [1]. The famous threshold theorems 
[2, 3], saying that arbitrary long quantum computation 
is possible provided the error per gate is below certain 
threshold has given the hope, that it is possible in prin- 
ciple to overcome the decoherence. However the initial 
theorems are based on phenomenological model of noise, 
and the problem, has not been solved within Hamilto- 
nian dynamics [4-7]. Even the problem of whether one 
can store qubits is open. 

There is, though a class of candidates for quantum 
memories, which are in between realistic description and 
the phenomenological one: the Kitaev models of topolog- 
ical quantum memory [8-10]. There is a heuristic reason- 
ing, according to which such memories are instable in two 
dimensions [8, 11], and stable in four dimensions (simi- 
larly like Ising model represent a stable classical memory 
in 2D, but not in ID) [8]. Behaviour of of Kitaev mod- 
els in finite temperature was then investigated (see e.g. 
[12-16]). Quite recently the thermal instability of 2D 
model has been rigorously proved in [17]. In the present 
paper, we deal with the 4D Kitaev's model of Ref. [8] 
and prove rigorously, within Markovian weak coupling 
approximation, that the model provides thermally sta- 
ble qubit. To this end we use the formalism of quantum 
semigroup theory [18], which has been successfully ap- 
plied to analysis of Kitaev 2D model in Ref. [17]. As a 
byproduct we obtain a very useful general upper bound 
for decay rate. We perform our analysis in parallel for 
3D and 4D case. Indeed, though in 3D case only one 
of the qubit observables is stable, as argued in [8], it is 
much more transparent and the reasoning is the same 
as in 4D case. Since the very stability of qubit is not 
sufficient for a good quantum memory, we also discuss 
the open problems concerning existence of self-correcting 



quantum memory. Implications for description of ther- 
modynamical limit are also discussed. 

The paper is organized as follows. In section II we 
provide some basic notions and results concerning Marko- 
vian weak coupling Hmit. We show, in particular, how 
the rate of decay expressed in terms of noise generator 
is related to fidelity criterion. Finally we provide a gen- 
eral upper bound for decay rate. In section III we show 
that analysis of noisy evolution of some particular topo- 
logical observables is reduced to the study of a classical 
model. Next (sec. IV) we provide conditions for stability 
of these observables in terms of one-step autocorrelation 
functions. In sec. V we finally prove the stability of the 
observables. In sec. VI we provide polynomial algorithm 
to measure the observables. Finally (sec. VII) we discuss 
remaining open problems for existence of self-correcting 
quantum memory, as well as importance of the result for 
description of systems in thermodynamical limit. 



II. MARKOVIAN APPROXIMATION IN WEAK 
COUPLING LIMIT 

Let us first we briefiy sketch the general setup and 
properties of Davies generators. A quantum system with 
discrete energy spectrum is coupled to a collection of heat 
baths leading to the global Hamiltonian 



H = H 



sys 



with 



(1) 

where the Sa are system operators and the /„ bath op- 
erators. The main ingredients are the Fourier transforms 
ha of the autocorrelation functions of the fa ■ The func- 
tion ha describes the rate at which the coupling is able 
to transfer an energy u from the bath to the system. Of- 
ten a minimal coupling to the bath is chosen, minimal in 
the sense that the interaction part of the Hamiltonian is 
as simple as possible but still addresses all energy levels 
of the system Hamiltonian in order to produce finally an 
ergodic reduced dynamics. The necessary and sufficient 
condition for ergodicity is [19, 20] 



(2) 
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i.e. no systom operator apart from the multiples of the 
identity commutes with all the Sa and H'^y^. 

We begin by introducing the Fourier decompositions of 
the Sa 's as they evolve in time under the system evolution 

e^*^"^ Sa e-'*^"^ =J2^- (a;) e"'-*. (3) 

Here the u are the Bohr frequencies of the system Hamil- 
tonian. From self-adjointness we have the relation 

Sahuj) = Sa{i0)l (4) 

The weak coupling limit procedure then returns the fol- 
lowing equation for the evolution of the spin system in 
Heisenberg picture 

^ = i[H^y^, X] + CaisiX) =: CiX) (5) 
^disiX) = i^^/i«(u;)(5iH[X,5„H]+ (6) 

a oj 

+ [SUuj),X]Sa{u;)) (7) 

For thermal baths one has moreover the relation 

ha{-uj) =e-^^'^ha{u;) (8) 

which is a consequence of the KMS condition [18]. The 
operator C generates a semigroup of completely positive 
identity preserving transformations of the spin system. 
It describes the reduced dynamics in the Markovian ap- 
proximation and enjoys the following properties 

• The canonical Gibbs state with density matrix 

is a stationary state for the semigroup, i.e. 

Tv(p/3e*^(X)) =Tr(p^X). (10) 

• The semigroup is relaxing, i.e. for any initial state 
p of the system 

^hm Tr(pe*^(X)) =Tr(p^X). (11) 

• Furthermore, the generator satisfies the detailed 
balance condition, often called reversibility. Writ- 
ing <5(X) := {n^y\X\, 

[5,£dis] = (12) 

and 

Tr(p;3yt£dis(X)) = Tr(p/3 (£dis(>^))^ x) ■ (13) 

The last equation expresses the self-adjointness of 
the generator with respect to the scalar product 



defined by the equilibrium state. The space of ob- 
servables equipped with the scalar product 

(X,y)^:=Trp^Xty (14) 

is called the Liouville space and the generator of 
the reduced dynamics is a normal matrix on that 
space, i.e. the Hcrmitian and skew-Hermitian parts 
of the generator commute. 

Finally it is known that — £ is a positive operator, hence 
it has nonnegative eigenvalues. Moreover i(7) = 0, and 
for ergodic systems eigenvalue is nondegenerate. 

A. Autocorrelation functions, decay rate and 
fidelity 

Suppose that for observable X satisfying (X, X)^ = 1 
(X, 7)^ = we have 

-(X,/:(X))^<e. (15) 

Then the autocorrelation function of the observable sat- 
isfies 

(X,e^*X)0>e-^* (16) 

One proves it easily, by decomposing X into normalized 
eigenvectors of £, and using convexity of function e~^. 
Thus to show that an observable X is stable, it is enough 
to estimate — (X, £(X)), which can be therefore called 
derM,y rate for the observable X. If this quantity decreases 
exponentially with size of the system, we obtain stability. 

Let us now rephrase it in the language of fidelity. 
Namely, suppose we have obscrvablcs X and Z satisfy 
commutation rules of Pauli algebra. They generate sub- 
algebra which defines a virtual qubit, the one to be pro- 
tected. Let the induced tensor product on the total 
Hilbert space be 

n = nQ® Jianc- (17) 

Now, we fix some state Panc on the system Hanc- For 
any state ijj of qubit the initial state of the total system 
is PQ,Qnc(0) = IV') ("01 ® Pane- Then the system evolves 
into state PQ,anc{i), and finally, the ancilla is traced out. 
Thus the fidelity is given by 

= {^\PQ\^m (18) 

where 

p°^'{f)=Tlanc{pQ.anc{t)). (19) 

Let us denote the fidelity averaged uniformly over the 
states of qubit by F. 

Proposition 1 With the above notation, suppose now 
that the Gibbs state is of the form 

P0 = ® Pane, (20) 
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where Panc is a state on ancilla. We then have 

F > e^'X)0 + {Z, e'-'Z)0) > e"^* (21) 

where e is upper bound for the rates —{X,jC{X))i3 and 
-{Z,C{Z))p. 

Proof. Let be given by 

F. = \{Fm)+F{m (22) 

where |0), |1) are eigenstates of X treated as observable 
on system Q. Similarly we define F^. Using results of 
[21] and [22] one finds that 

F>F^ + F,-l (23) 

Thus it is enough to estimate e.g. F^- Using the property 
(20) and orthogonality X A- I one finds that 

F, = i(l + (X,e^*X)^). (24) 
Combining the last two formulas with (16) ends the proof. 

B. Upper bound for decay rate 

We now present a useful bound for decay rate, which 
holds for operators X which are eigenvectors of [H, ■] . 
For such operators, one computes 

- {X,C{X))0 = Y,K'^){{[Sa{uj),X],[So,{uj),X])0 

+ e-^''{[S^{~uj),X],[S„{-Lo),X])f^ < (25) 
<2Y,H^){{Sa{co),X],[Sc.{iv),X])0 < 

< 2Va,^([5a(a;),X],[5„(a;),X])/3. (26) 

where 

ftmax = sup/i(u;). (27) 

Since X and Sa{Lu) arc eigenvectors of [H,-] it follows 
that [Sa{ijj),X] are eigenvectors of [H, ■] too, hence they 
are mutually orthogonal. We thus can write 

J2{[Sc {iv), X] , [5„ (u;) , X])^ = J2{{Sa {io) , X] , [S^ {cj'), X]) 

(28) 

However from definition of Sa{oi>) it follows that 

Y,SU^)=S^ (29) 

This gives 

-{X,£{X))0 < 2hn.^Y^{[Sc,,X],[Sa,X])f3 (30) 

a 

The advantage of the formula is that the only place where 
the self-Hamiltonian appears is the Gibbs state in scalar 
product. 



III. FROM QUANTUM TO CLASSICAL IN 
KITAEV-TYPE MODELS 

We consider a system of A'' spin-1/2 systems. For any 
set S of spins let us denote Xs = Hj^s^'j , Zs = Hj^so'j- 
Consider now Hamiltonian of the form 

H=-J2x,-J2Zc (31) 

s c 

and we assume that the sets s and the sets c are chosen 
in such a way that the operators Xg and Z^ commute 
with each other. Consider also the following couphng to 
environment 

Hint = 1] CT," ^ fj + E ® fj ■ (32) 

3 3 

Then Davies operators fall into two types: 

fla = CliPc. (33) 

ba = criRa (34) 

where Pa belong to algebra spanned by those operators 
Zc whose support contains the j-th spin and Ra belongs 
to algebra spanned by operators Xg , whose support con- 
tains j-th spin. (If the spin does not belong to support 
of any s, then Pa = I, and similarly for R. However, in 
Kitaev-type models this latter case does not occur). The 
dissipative generator has the following form 

L = C^+C„ (35) 

where Cx, Hz consist of Davies operators of type a and h 
respectively. The Davies operators describe the elemen- 
tary noise processes. In 2D model, they are creation, 
anihilation and motion of two types of point-like anyons. 
In 4D model, excitations are not point-like, and the pro- 
cesses are creation, anihilation and two types of modifi- 
cation of loops (see [8], sec. X, and sees. V, III A of the 
present paper). 

Consider now observables of the form Xs and Zt, 
where 5, T are some subsets of spins. Let us assume 
that Xs and Zt commute with all X^ and Z^. Then Xs 
commutes with Davies operators of type a and Zt com- 
mutes with Davies operators of type b. Therefore from 
(7) we get that 

£(Xs)=£.(Xs), £(Zt) = £.(Zt) (36) 

Consider now a modification of the model. Let the 
Hamiltonian be of the form 

H=-Y,X, (37) 

S 

and the coupling with environment be of the form 

H,nt = Y^(yi®h- (38) 

3 

Then dissipative generator for this model consists of 
Davies operators (34) i.e. it is given just by C^. We 
obtain 
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Proposition 2 Let Xs = Ilj,=s(7j , Zc = TljecO'j where 
the sets s and c are chosen in such way that Xg and 
commute for all s,c. Consider Xs which commutes with 
all Xg and Z^. Then 

C{Xs) = C'{Xs) (39) 

where C is dissipative generator coming from 

H = -Y,Xs-Y.^c, hm = Y.(t^<s) fj + ^ a| ® . 

s c j j 

(40) 

and CJ is dissipative generator coming from 

H'=-Y,Xs. m^, = ^a]^f,. (41) 

s j 

Moreover 

TV(p^Xt£(X)) = Tv{p'pX^C'{X)) , (42) 

where Pjs = ^e~^^ and p'^ = -^e^^^ respectively. 
Analogous result holds for Zt, which commutes with all 
Xg and Zc- 

Remark. Further in text, {■,■)() will denote scalar 
product with the Gibbs state of type (with suitable 
H', depending whether we talk about X or Z). 

A. Observables X and Z 

1. 3D Kitaev model 
The Hamiltonian for 3D Kitaev model is given by [8] 
H = -J2Xs-J2^'^ (43) 

s c 

where each s denotes set of four plaquettes which share 
common Hnk, and and each c is six plaquettes forming 
cube. We will now define a class of observables of interest. 
To this end we will use observable Xc with C being set of 
parallel plaquettes forming a loop that winds around the 
torus (there are three homologically inequivalent choices, 
we will consider a fixed one of them) . Such observable is 
very unstable, hence we may call it "bare qubit observ- 
able". One needs to "dress" it with another dichotomic 
observable which would store the error syndrome. The 
latter observable will then belong to the abelian algebra 
spanned by star observables Xg, hence depending solely 
on atomic projectors of the algebra which correspond to 
configurations JC of excited finks (stars can be labeled by 
the links - their centers). Let us call the projectors P/c- 
The needed observable will be thus of the form 

F.=Y,^'cPjc, (44) 

K 

where \)c = ±1. We shall not determine the values of 
Ax; at the moment. They will emerge from our analysis 



of stability in sec. V and will be then described in sec. 
VI. 

The full "dressed observable" is the product XcF^. 
According to Proposition 2 it evolves according to clas- 
sical model with Hamiltonian 

Hx = -J2^- (45) 

coupled to environment via operators (t|. The model is 
known as Z2 gauge Ising model (the Ising variables are 
in our case eigenvectors of trj) [23]. 

One can define analogous observable ZpF^. However 
in 3D there will be no X-Z symmetry. The observable Zp 
is associated with plane, and atomic projector of algebra 
spanned by Z,. is labeled by configurations of points (i.e. 
the plaquettes) rather than by loops. Observable ZpFz 
is evolving according to the model with Hz = — X^c 
coupled via crj. It wifi not be stable (as pointed out in 
[8]) and most likely, one can prove it by use of techniques 
worked out in [17]. 

2. 4:D Kitaev model 

In four dimensional model the spins again sit on pla- 
quettes, and the Hamiltonian is similar as in 3D case: 

H=-J2Xs-J2Zc (46) 

s c 

The only difference is that the star s has six plaquettes, 
because there is six plaquettes common to a single fink. 
Thanks to it there is symmetry: We fix two planes pi and 
P2 on the lattice and on the dual lattice, respectively, ob- 
taining bare qubit observables Xp^ and Zp^. Then can- 
didates for stable observables will be the dressed ones 
Xp^Fx, Zp^Fz. The latter will again evolve separately, 
and since 4D lattice is self-dual, the evolutions are the 
same. We arrive at the 4Z) Z2 gauge Ising model. 

If we prove that e.g. observable of the form Xp^F^ is 
stable, then also similar Zp^F^ will be stable too, so that 
we wifi obtain stable qubit. 

IV. STABILITY CONDITIONS FOR KITAEV 
MODEL 

A. Bound for decay rate for dressed observables 

The bound (30) applied to generator consisting of 
Davies generators (33), (34) takes the form 

-{A,C{A))0 < 2VaxI]([a;, A], [a],A])0 + 
j 

+ ^([a|,A],[<7|,A])^. (47) 
3 

The quantity /)-max given by (27) is a constant indepen- 
dent of the size of the system. This is due to the fact that 
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Kitaev models exhibits strong locality property, implying 
that there is a constant number of frequencies involved 
in the generator (e.g. just one positive frequency in 2D 
model) which are independent of the number of spins N. 

Since the observables Z — ZpFz, X = XcF^ (or anal- 
ogous ones from 4D model) commute with Hamilotnian, 
the bound is applicable. We obtain 

-{X,C{X))p < 4/i^ax^(l - 

3 

-{Z,C{Z))p < 4h^,^Y^{l - {Z,a]Za])0) (48) 
j 

where j runs over all spins. We see that the problem of 
decay of time autocorrelation function has been reduced 
to the much simpler problem of "one step" autocorrela- 
tion function. 



B. Gibbs state is concentrated on configurations 
without long loops 



where i5 = /3 — In /x is positive and does not depend on 
the size of the system. We then evaluate probability of 
appearing a configuration that has a loop greater than 
L' 

P{1 > L') < poly{L) e-'' = poly{L)e-''^' 

l=L' ^ 

(53) 

Thus we see that below Tcrit the probability of obtain- 
ing e.g. a loop of length L/8 or greater is exponentially 
decaying in L. 



C. Stability of Kitaev 4D model 

In next section we shall prove that for configurations 
having only loops shorter than L' = L/8 a single flip does 
not change observables X and Z for Kitaev 4D model. 
This implies that 



First wc will estimate probability that a configuration 
has loop of length I. We shall use the Peierls argument 
following Dennis et al. [8] and Griffiths [24]. To this 
end we first estimate probability that a fixed loop A with 
length I emerges. Let C be the set of all configurations 
which contain loop A. The probability is then given by 



P(A) = 



(49) 



where in denominator we have sum over all configura- 
tions. For any configuration /C containing A we fiip spins 
on a chosen surface whose boundary is A, obtaining new 
configuration /C* which differs from /C only in that the 
loop A is not present anymore. Hence i?(/C) = E{JC*)e~^^ 
(or the quantities here are taken to be dimensionless). 
Thus we write 



P(A) = 



-/3/y- -I3E(K.') 



-I3E{K) 



(50) 



Leaving in denominator only configurations /C*, we can 
only decrease it, so that P{X) < e^^^ . 

Now, the probability P{1) of appearing a configuration 
which has a loop of length I is bounded by the number 
of all possible loops of length I times e"'''. A trivial 
bound for the number of loops in cube of linear size L in 
dimension d, that start from a fixed node is 2d{2d — 1)'. 
This should be multiplied by the number of nodes, which 
is proportional to the volume i.e. polynomial in linear 
size L of the system. Finally, wc obtain that 



P{1) < poly{L)ii^e 



-pi 



■ poly{L)e 



-l([i-\nt,) 



(51) 



where is a constant depending only on d. Thus below 
certain critical temperature Tcrit we have 



P{1) < poly{L)e- 



(52) 



^(1 - {Z, a^Za^)p) < ^ 2P(Z > L') (54) 



so that 



-{Z,C{Z)) < poly {L)e- 



-5'L 



(55) 



where 5' = (5/8 is a constant that is positive below some 
critical temperature. The same happens for observable 
X, hence due to proposition 1 the decay time of fidelity 
is exponentially long in size of the system. 



STABILITY OF TOPOLOGICAL 
OBSERVABLES 



In previous section we have shown that below certain 
critical temperature T^it the Gibbs state is concentrated 
on configurations with short loops. Thus if on such con- 
figurations an observable does not change under single 
spin fiip, it is stable within the classical model. If in ad- 
dition it is of the special form XcF^ , then it is also stable 
within the quantum model (see section VC). 

In this section we shall build such observable. To this 
end we shall first define homology classes of spin con- 
figurations corresponding to configurations of loops with 
short loops only. We will then show that, as expected, 
single spin fiip does not change those homology classes. 
This impHes that any observable which depends solely 
on the homology classes docs not change under single 
spin fiips (for configurations containing only short loops) . 
This result holds for torus of any dimension. Subse- 
quently we shall show, that some observables of the form 
XcFx share this property. 
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A. Observables depending only on homology 
classes 

Let us introduce some notation. By S we will denote 
configuration of spins on the lattice (in tlic form of con- 
figurations of bits whose values encode spin orientation) . 
Given two spin configurations Si and 5*2, we can add 
them to obtain new configuration S. We denote it by 
S = Si (B S2, and the addition is bit-wise, modulo 2. 
I.e. if at given site the spins are the same, resulting spin 
is down, if they are different, resulting spin is up. We 
denote by 5*0 configuration of all spins down. 

By /C we will denote set of excited links. A link is 
excited, if the parity of spins on adjacent plaquettes is 
odd (in 3D a hnk has four such plaquettes, and in 4D - 
six ones). One finds that /C is sum of disjoint loops Ij 
(the loops can have self crossing at nodes) : 

Lemma 1 For given configuration K consider a con- 
nected set of links. It is sum of dosed loops, which visit 
each link and each node at most one time. Equivalently, 
it is a closed walk, which visit one link only at most once. 

Proof. The proof is by induction. 

We will call such connected sets "loops". We will say 
that a loop is short, when its length is no greater than 
cL, where c is a fixed constant, which we can take e.g. 
1/8. 

Any spin configuration S defines link configuration /C. 
We will then write S{1C). Of course for given /C there are 
many spin configurations leading to them. Sometimes for 
given S the corresponding K, will be denoted by dS and 
called boundary of S. 

Definition 1 By continuous deformation of spin config- 
uration we mean operation, which can he composed of the 
following elementary operations: flipping spins on all pla- 
quettes belonging to an elementary d-dimensional cubes. 

Remark 1. Continuous deformation does not change 
the configuration of links. For 3D easy to see: indeed, 
fiipping spins on faces of cube, change at the same time 
spins on two plaquettes adjacent to a link from the cube. 

Definition 2 We say that Si and S2 with empty bound- 
ary are homologically equivalent if they can be trans- 
formed into one another by continuous deformation. S 
is called homologically trivial, if it can be continuously 
transformed into Sq. 

Definition 3 We say that Si and S2 which have the 
same boundary are homologically equivalent and denote 
it by Si ~ 52, if Si ® ^2 is homolgically trivial 

Definition 4 ("Shortest configuration") Consider given 

IC = Uj • For each loop Ij fix a shortest surface whose 
boundary is Ij. Consider then S*{lj) which has spins up 
on this surface and all other spins down. The configura- 
tion S* = (BjS*{lj) will be called shortest configuration 
for C. 



Fact 1 All shortest configurations S* for given IC are 
homologically equivalent, provided IC contains only short 
loops. 

Proof. Take two different shortest configurations. We 
have 

s*i®s; = (s,[s*iii,)®s;{i,)] (56) 

However, each configuration 5'i(^j)©S'|(/j) is trivial. In- 
deed, since loop Ij is short then \Si{lj)\ and |«S'2(/j)| are 
small, and cannot form homologically nontrivial surface. 
I 

Definition 5 For JC containing only short loops, with 
any S leading to K. we can associate the homology class 
ofS®S*{K:). Denote it by h{S). 

Remark 2. For fixed IC obviously Si ^ S2 iff h{Si) = 
h{S2). Thus the above definition allows to ascribe labels 
to homology classes of spin configuration, by relating to 
distinguished class i.e. the class of S*. But the homology 
classes are now defined for any /C. Thus we will be able 
to ask later, whether a spin flip (which of course changes 
IC) can preserve homology class. For any IC there are 
eight homology classes in 3D case, associated with three 
possible ways of winding around torus. In 4£) there is 16 
classes. 

We have obvious fact: 

Fact 2 We have Si ® S2 = crjiSi) ® aj{S2), where aj 
flips j-th spin. 

Now we will show that for short loops, single spin flip 
does not change homology class of 5. To this end we first 
prove the following lemma 

Lemma 2 For IC containing only short loops we have 

a,{S*{IC))^S*{a,{IC)). (57) 

Here (Ji{IC) is understood as the configuration of loops 
arising from configuration IC by applying Uj 

Proof. Divide K, into two sets: Ki consisting of loops 
that contain some links from i-th plaquette, and IC2 
which does not contain Hnks from this plaquette. Then 
(Ti(/C) = CTiiJC) U ICi hence 

cJi{S*{IC)) = ai{S*{ICi)) e S*{IC2) (58) 

and 

S*{ai{IC)) = S*{ai{ICi)) 5*(/C2). (59) 

Thus only /Ci is in the game: 

ai{S*{IC))®S*{ai{IC))=ai{S*{ICi))®S*{ai{ICi)) (60) 

and therefore we have to show that right-hand-side of 
the above formula is homologically trivial. Indeed, the 
set ICi contains at most two loops independently of di- 
mension. Now, since loops are short, both ai{S*{ICi)) 
and <S'*((Ti(/Ci)) are small, and added together must give 
a trivial surface. | 

Now we are in position to prove the main result of this 
section 
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Proposition 3 Consider configuration of spins S for 
which IC has short loops only. Then single spin flip does 
not change the homology class of S. More explicitly, we 
have 

ai{S)®S*{ai{IC)) ^S®S*{K:) (61) 

Proof. By lemma 2 we have 

ai{S*{IC)) ~ 5*(a,(/C)). (62) 

By fact 2 we have 

ai{S)®ai{S*{IC) = S®S*{K:) (63) 

Combining the above two equations, we obtain the 
claim.| 

Thus any observable T which for configurations IC con- 
taining only short loops depends only on homology class, 
i.e. 

T{S)=T{h) (64) 

is dynamically stable within the model of Proposition 2 
below some critical temperature. 

B. Construction of stable topological observables 

Our bare observable will be Xt — IljgTf^ where T is 
chosen in such a way that Xt is invariant under flipping 
spins on plaquettes from any cube (i.e. it is invariant un- 
der continuous transformations). Examples of such ob- 
servables exists, as will be shown later in next subsection. 
We will show that one can find dichotomic observable E,. 
which will depend on given configuration S only through 
IC, such that the dressed observable XtF^ depends only 
on homology of S (for short loops) i.e. it is stable within 
the classical model. 

We begin with the following lemma 

Lemma 3 The observable Xt which is invariant under 

flipping spins on plaquettes of any cube is constant on 
homology classes for any fixed link configuration IC con- 
taining only short loops (cf. definition 5). 

Remark. Note that this does not mean that Xt is sta- 
ble. Indeed, for any fixed link configuration, it is con- 
stant on the whole homology classes. However if the link 
configuration changes, it may change sign on the same 
homology class. The stable observable described in pre- 
vious subsection has the same value on a given homology 
class independently of link configurations, provided there 
are only short loops. 

Proof. Consider arbitrary spin configurations S and 
S' whose boundary is IC, and which are in the same class 
of homology, i.e. Si = S ® S* \s homologically equivalent 
to S2 = S' ® S* . Therefore Si can be transformed into 
^2 by flipping spins on a set of elementary cubes. This 
does not change the sign of Xt, so that Xt has the same 



sign on Si and S2. Thus it has the same sign also on S 
and 5". | 

Now we are in position to build stable observable. Now 
let us assume that 

Xt{Si)Xt{S2) = Xt{Si e 52). (65) 

We stress here that this assumption is easily seen to hold 
for particular observables considered in next subsection. 
(One can actually show, that it is true in general for 
observables satisfying assumptions of the above lemma) . 
Using this we can write 

Xt{S) = Xt{S ®S*® S*) = XTiS © S*)XTiS*) (66) 

Since homology class of 5"* is always the same for short 
loops (independently on possible amiguity of S* for given 
loop) then Xt{S*) depends only on the loops configu- 
ration: Xt{S*) = X"{IC), so that Xt{S) = Xt{S © 
S*)X" (IC). Now, since for fixed loops configuration Xt 
depends only on homology class and the loops configura- 
tion for S (B S* is always null (as S (B S* does not have 
a boundary), we get that Xt{S (B S*) depends only on 
homology class of S (B S* . Therefore, according to def- 
inition 5, it depends only on homology h of S. Hence 
Xt{S © S*) = X'{h) and we have 

Xt{S) = X'ih)X"{IC). (67) 

Then the following observable 

T{S) = Xt{S)X"{IC) (68) 

depends only on h. The above observable is defined un- 
ambiguously only for spin configurations leading to short 
loops configurations. This is because X" is only well de- 
fined only on short loops configurations. We then extend 
the definition of T to all spin configurations, by letting 
X"{IC) = 1 for all other loops configurations. Thus we 
shall take = X" and obtain that XtF^ depends only 
on homology of spin configuration, hence is stable within 
classical model. 



C. Observable stable within quantum model 

The observable constructed in the previous subsection 
is stable within classical model, because it depends only 
on homology class. However, we know that only special 
observables from the classical model evolve in the same 
way in quantum model. E.g. the observables of the form 
XtF^, where F^ is from algebra generated by star op- 
erators Xs,& and Xt commutes with Zc. Here wc shall 
focus on construction of Xt since it determines F^ via 
considerations of the previous section. 

Now, let us note that the first condition means simply 
that Fx depends only on loops. The second condition 
means that Xt does not change under fiips on all pla- 
quettes of an elementary cube. Thus the observable (68) 
is of the above form, hence it evolves in the same way 



both in quantum and classical model, hence it is stable 
also within quantum model. 

The last thing is to assure that the observable T is non- 
trivial, i.e it is not identity. To this end we have choose 
the set T in a special way, such that on spin configura- 
tions without boundary, Xt can take different sign. 

For 3D it will be nontrivial loop in dual lattice, i.e. 
straight Hue consisting of parallel plaquettes. The fact 
that it is loop in dual lattice, implies that it Xt is in- 
variant under continuous transformations. Since it is 
nontrivial, then Xt have value —1 for spin configura- 
tion consisting of plane of fiipped spins perpendicular to 
T, while it takes value 1 on homologically trivial spin 
configurations. Since there are three possible choices of 
inequivalent nontrivial loops, we can construct three in- 
dependent observables. 

In AD we take T to be plane in dual lattice, i.e. the 
value observable Xt is defined as a product of values of 
all plaquettes belonging to the plane T. Again, Xt does 
not change under flipping spins on cube because arbitrary 
cube has exactly two plaquettes in common with such a 
plane. For this reason it will be also 1 on homologically 
trivial spin conflgurations. However it will take value —1 
on the conflguration consisting of flipped spins on a plane 
T' whose intersection with T is a single plaquette. Note 
that since there are six homologically nontrivial planes, 
we can construct six independent observables of this sort. 

Now, since the torus in AD is selfdual, we can consider 
dual observable i.e. % = Z^F^, and depends only 
on conflguration of three dimensional cubes (such cubes 
are dual to link). Since Fz and F^ commute, and planes 
T and T' intersect only in a single plaquette, we obtain 
that Tz and Tj anticommute, so that they form a qubit. 



VI. POLYNOMIAL ALGORITHM FOR 
MEASURING THE TOPOLOGICAL 
OBSERVABLES 

The observables are symmetric, so it is enough to show 
algorithm for one of them, say The algorithm is the 
following. 

1. Measure all spins. 

2. Multiply outcomes on a flxed plane in dual lattice, 
this gives "raw value" of the observable. 

3. Identify the loops. 

4. For "short" loops we identify associated surfaces 
(the ones homologically equivalent to shortest 
ones). 

5. If an odd number of surfaces crosses a flxed plane 
in dual lattice, multiply the "raw value" with —1. 

The step 2 corresponds to measuring the bare observ- 
able Xt, while the steps 3-5 deflne observable F^. The 
multiplication in last step produces the stable, dressed 
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FIG. 1: Efficient algorithm for determining surface closing the 
loop. 

observable Tx = XtFx- The only nontrivial problem 
here is to argue that the step 4 is polynomial. It is ac- 
tually enough to show that for a flxed loop, one can flnd 
efflciently a surface which is contained in the smallest 
cube containing the loop. 

Moreover, it is enough to flnd a protocol which in ef- 
flcient way allows to flnd spins which, if flipped, reduce 
length of the loop by some amount (in our protocol, it 
will be reduced by two). 

The protocol is the following. We flrst choose a Carte- 
sian frame. We start with a link of the loop, and move 
along the loop. If there is ambiguity (the loop crosses it- 
self) the priority is set by the chosen frame: if only we can 
we go in positive direction of the axis with the smallest 
number. If not, then we go in negative direction of the 
axis with the smallest number. The same rule governs 
choice of the starting link. 

The walk is stopped, if we are forced at some point to 
go in opposite direction to any of the previous steps (see 
flgure) . 

When the walk is stopped, the link at which we stopped 
and the last "opposite" link, determine uniquely the set 
of plaquettes. This is because all the links of the walk 
lying between two "opposite Hnks" are perpendicular to 
them. Now, after flipping spins on the set of plaquettes, 
the two opposite links are removed from the curve. Note 
that this flipping may further diminish the length of loop, 
if by a chance, the chose plaquettes have some other Hnks 
common with the loop. It may also divide the loop into 
smaller ones, however their joint length is not longer than 
I - 2. 



VII. CONCLUDING REMARKS 

We have shown that within Markovian weak coupling 
approximation, there exist a stable quantum subsystem 
in four dimensional Kitaev model of [8]. While the qubit 
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is indeed stable, there are several other drawbacks, which 
makes the question of existence of self correcting quan- 
tum memory still open. Minimal requirement for good 
quantum memorj' is that it should allow to encode arbi- 
trary state of qubit (encoding) , then to store it for long 
time (storage) and finally perform a measurement in arbi- 
trary basis (readout). It would be also good if the mea- 
surement is repeatable. The present result shows that 
storage is possible, but docs not touch the problem of 
preparation and measurement. Actually, the algorithm 
for measuring topological observables is highly destruc- 
tive, hence non-rcpcatablc. The encoding and read-out 
one usually performs by preparing qubit in a standard 
state, and also measure standard observable, the rest be- 
ing done by gates. Also repeatability can be then assured, 
if one can perform c-not gates on tlic protected qubits. 
However the problem with the Kitaev's model is that it 
does not support universal computation. A possible so- 
lution of this problem is to use the version of topological 
quantum memory developed by Bombin and Delgado [25] 
which supports universal computation (we shall present 
the dynamical analysis of these models elsewhere) . How- 



ever, still there is a separate problem of preparation of 
the qubit in standard state. 

Finally, let us mention, that the present result has sep- 
arate implications in statistical physics. In the standard 
approach to large quantum systems the metastable states 
of encoded qubits like those found for Kitaev models dis- 
appear in the thermodynamic limit merging into a clas- 
sical simplex of equilibrium (KMS) states [11]. On the 
other hand they carry an interesting topological struc- 
ture which might be physically relevant. In this context 
it is interesting to ask for a new description of infinite sys- 
tem, which would take into account such new metastable 
states. In particular, phase transitions which lead to such 
curious states require further investigations. 
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